% fig2b.m ——— 最大李指数谱 ———
clear; clc; close all;

%% —— 参数设置 —— 
beta   = 1;    xi    = 0.8;    mu    = 1;
A      = 3.4;  omega = 1;

% —— 时间控制参数 —— 
dt    = 0.01;    % 时间步长
T0    = 200;    % 丢弃瞬态时间
Tsim  = 400;    % 统计时间

% —— 扫描 α —— 
alpha_min  = 1.0;
alpha_max  = 4.0;
Nalpha     = 200;
alpha_list = linspace(alpha_min, alpha_max, Nalpha);

LEs_list = zeros(1, Nalpha);

%% —— 并行计算最大李指数 —— 
parfor ia = 1 : Nalpha
    alpha = alpha_list(ia);
    % 将 dt、T0、Tsim 传入
    Ls = LEs(alpha, beta, xi, mu, A, omega, dt, T0, Tsim);
    Ls_list(ia) = Ls;
end

%% —— 绘图 —— 
figure;
plot(alpha_list, Ls_list, 'b-', 'LineWidth',1.5);
hold on;
plot(xlim, [0 0], 'k--', 'LineWidth',1);  % LE=0 参考线
xlabel('\alpha','FontSize',12);
ylabel('LE 指数','FontSize',12);
title('Fig.2(b) LE 指数图','FontSize',14);
xlim([alpha_min alpha_max]);
grid on;
